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After inflation, a period of preheating may have produced a stochastic background of high fre¬ 
quency gravitational waves (GWs) that would persist until today. The nature of the inflaton’s cou¬ 
pling to Standard Model or other fields is unknown, so it is useful to ask what features such fields 
may typically have, and how these affect predictions for the GW’s produced. Here we consider the 
inflaton to be coupled to a light scalar field, and show that even a very small quartic self-interaction 
term will reduce the amplitude of the GW spectrum. For self-coupling A x > g 2 , where g 2 is the 
inflaton-scalar coupling, the peak energy density goes as f2gw x V^gW x_ °^ ~ (<7 2 /A x ) 2 . A conse¬ 
quence is that if the universe reheats through an inflaton-Higgs coupling then the spectrum would 
be suppressed but the dynamics would be sensitive to the Higgs potential near the energy scale of 
inflation. 


I. INTRODUCTION 

Inflation leaves the universe cold and nearly empty of particles, so there needs to be a reheating mech¬ 
anism for energy transfer between inflaton and Standard Model fields in order to create the thermalized 
particles that existed before Big Bang Nucleosynthesis began. This is typically modeled by a small, direct 
coupling between inflaton and another field. The first discussions of reheating [1-7] studied a perturbative 
calculation of inflaton decay into the coupled field, with energy gradually transferred to matter fields. (Also 
see the reviews [8, 9].) 

However, inflaton decay occurs in the context of large, coherent field oscillations and nonperturbative 
effects should also be taken into account [6, 7, 10, 11], Typically, the inflaton cf) is considered to be coupled 
to a field x by an interaction \g 2 4> 2 x 2 , which is y ’s only potential energy term. As the inflaton oscillates 
about the bottom of its potential after inflation, the phenomenon of parametric resonance leads to some 
modes of the decay product \ being excited at an exponential rate. This effect, which may occur briefly at 
the beginning of a longer period of reheating, is called preheating. (Most of the work on this subject has 
been in the context of direct couplings between inflaton and matter fields; see [12] for a scenario that does 
not require this.) 
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Preheating in these models can produce gravitational waves [13-18], since the exponential amplifica¬ 
tion of certain modes leads to a large contribution to anisotropic stress, which sources tensor perturbations. 
Predictions for the resulting spectrum are around h? fl gw ~ 10“ 10 and / ~ 10 4 to 10 6 Hz today for massive 
or A0 4 inflation or could be as low as 10 2 to 10 3 Hz for hybrid inflation models. Some work [17, 19-24] 
has addressed this problem in the context of various models that relate to processes that are more specific. 
These find a wide range of possibilities. For example, [21] found that decay into fermions after inflation 
could produce Q gw ~ 10 -12 to 10” 18 , / ~ 10 9 to 10 10 Hz today, depending on the parameters. 

These tend to fall outside the range of current, planned or proposed gravitational wave experiments 
such as Advanced LIGO and VIRGO, KAGRA, Einstein Telescope, eLISA, DECIGO or BBO (for an 
exception, see [20]). Roughly speaking, these are most sensitive to frequencies around 10~ 3 to 10 3 Hz 
and signal strength corresponding to // 2 O gw ~ 10 5 to possibly 10 14 . (See [25, 26] or the review [27]. *) 
LIGO and VIRGO have jointly placed upper bounds on a stochastic gravitational wave background on the 
order of O gw ~ 5 x 10 -6 around 10 2 Hz [28]. Gravitational wave detection at MHz frequencies has also 
been considered [29-31]. It has not been a major focus, though, since comparatively reliable astrophysical 
sources (e.g. neutron star mergers) are not expected in this frequency range. 

This motivates the study of how robust are the predictions for the gravitational wave spectrum from 
preheating. We would expect that a realistic preheating process in the early universe would include cou¬ 
plings of the decay product to other fields, as well as possible self-interactions. It will be useful to know 
whether these can significantly affect the observability of such a process. 2 Specifically, it would be inter¬ 
esting to answer the question “Given a model of preheating with some self-interaction strength, how does 
one estimate the overall gravitational wave production?” This is analogous to the discussion in [33], which 
estimates the maximum energy density in gravitational waves that could be produced by a cosmological 
process such as preheating. 

Previous work has shown that for self-couplings A x ~ C9(10 -2 ) V> g 2 , where g 2 is the coupling 
between the inflaton and scalar, parametric resonance can be significantly affected [34, 35]. However, there 
has been little discussion of gravitational wave production in this scenario. 3 Therefore it is difficult to give 
a thorough answer to the above question based on the existing literature. This also means that it is unclear 
how general a gravitational wave prediction is when it ignores interactions of the decay products. 

In this work, we begin to address this by studying the development and termination of parametric 

1 Note that some results are given in terms of /i 2 f2 gw , others in terms of fi gw and still others in terms of strain h, which is 
distinct from today’s Hubble constant in units of 100 km/s/Mpc that appears in h 2 fl gw . Consistent comparison of experimental 
sensitivities is discussed in [25], 

2 While this paper was in preparation, another work [32] appeared that addresses some of these questions. We will discuss it in 
Sec. VI. 

3 A study of gravitational waves in M-flation preheating [36] mentions that a self-interaction can suppress the resonance, but does 
not quantify this in a way that allows comparison with [34]. 
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resonance and the production of gravitational waves in the context of Ad> 4 chaotic inflation coupled to a 
self-interacting light scalar field. We verify by lattice simulation that the resonance terminates early for 
self-coupling A x > g 2 , demonstrating the condition p® nal ~ /y 2 / A ;v mentioned in footnote 19 of [34] (their 
g is our g 2 ), and show that this leads to significant suppression of gravitational wave production. The res¬ 
onance terminates early because the self-interaction term allows more efficient rescattering of particles out 
of the resonant mode, and this can be characterized by a condition comparing the energy density associated 
with the self-interaction to the inflaton-scalar interaction energy. The early termination of the resonance 
means that there is less energy in the light scalar’s fluctuations, which directly source gravitational waves. 
Therefore, gravitational wave production is reduced. For A x > A * = g 2 , the energy density goes as 

^gw x) ~ (fi , 2 /A x ) 2 ^gw x=0) - 

In Sec. V we show that this result is robust to changes in initial conditions, and that the same scal¬ 
ing occurs in massive (m 2 c/) 2 ) inflation. Although this suggests generality to inflationary models that are 
quadratic or quartic about the minimum, we point out that an important goal of future work is to understand 
the effect of realistic interactions on other models that have predicted gravitational wave spectra. 

As an application of this result, one could imagine the universe reheating by a coupling between the 
Higgs and inflaton, and we argue in Sec. VI that such a scenario would likely produce no observable grav¬ 
itational radiation. This is due to the size of the Higgs self-coupling, despite its eventual running to zero 
in the Standard Model. However, we point out that even a resonance too brief to produce observable grav¬ 
itational waves could be relevant for the issue of vacuum stability. Finally, if the inflaton preheats a scalar 
field with an extremely small self-coupling, then the gravitational wave spectrum could directly measure 
this potential. 


II. MODEL 

Representing the universe by a spatially flat Friedmann-Robertson-Walker metric, we will describe 
gravitational waves as transverse and traceless perturbations to this metric, specifically as hij such that 

ds 2 = a 2 (g) (—dr/ 2 + (Sij + hij)dx l dx^) (1) 

with dihij = 0 and ha = 0. We will take the inflaton to be a real scalar field, <p(t, x), and consider it to be 
coupled to a massless real scalar field \(f, x), with potential given by 

V = ^Xcj) 4 + ^A x x 4 + (2) 

Here we have chosen to study A0 4 chaotic inflation, and this requires some justification since standard slow- 
roll inflation with this potential is inconsistent with Cosmic Microwave Background (CMB) observations 
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[37]. Much literature on gravitational waves from preheating takes the potential as jA</> 4 , in particular the 
thorough numerical study [16], whose model corresponds to ours with the choice A x = 0. We expect 
the qualitative nature of our results to be relevant to a broad range of inflationary scenarios (this will be 
discussed further in Sec. VI), and it will be useful to refer to specific previous results in order to understand 
the production of gravitational waves. 

We are also studying the behavior of a “light” scalar field, and so we neglect a x mass term in com¬ 
parison with the effective x mass that comes from the interaction term kW- Comparing these terms 
using the amplitude of the 4> oscillations shows that this is roughly equivalent to requiring the x mass to be 
rn x <C g 2 / A x 10 12 GeV. 

Here the inflaton self-coupling is set by the amplitude of the scalar power spectrum of the CMB as 
A = 10 -13 . The unknown coupling g 2 must be small, but we will also take it to be larger than A; in 
terms of the resonance parameter q = g 2 /A this means 1 « q <C A -1 ; here we will examine the range 
10 Q < 2000, which contains most of the region with the largest gravitational wave production [16]. We 
will see that this peaks around q ~ 100 — 200 and falls off slightly as q gets larger or smaller (see Fig. 2d), 
although there are examples with smaller q that do not exactly follow this trend [16]. We consider the light 
scalar’s self-interaction in the range A < A x < 1. 

We study the dynamics in this model beginning at the end of inflation, to = 0, once the comoving 
horizon (aiF) -1 begins to expand, with the inflaton as a homogeneous field given everywhere by <po = 
0.342 Mpi. 4 The field y is a light “spectator” field during inflation, and at the end of inflation each y mode 
is in the de Sitter vacuum state. As shown in previous work [38, 39], as the inflaton decays the quantum 
state quickly approaches a semiclassical regime with large occupation numbers, and the evolution here is 
equivalent to the classical evolution of an initial classical distribution that gives 

(|Xfc(0)| 2 ) = 1/(2A 3 / 2 0oCjfc), Xk( 0) = (iwfc + H( 0)) y fc (0) (3) 

at the beginning of reheating. 5 The dynamics considered here occurs on sub-horizon scales. 6 

Since cf) is homogeneous, the equations of motion for these fields in a spatially flat Friedmann- 

4 This particular point along the inflaton’s phase space trajectory is identical to that of [16]. This choice is further addressed in 
Sec. V. 

5 Xk and cvk are defined below. The specific implementation for initial field conditions of [38] is as described in the documentation 

for LATTICEEASY code, available at http://www.felderbooks.com/latticeeasy/. 

6 For the typical example q = 120, numerical results show that preheating begins at about H = 1.1 x 10 -9 Mpi and a = 5.5 

(for a = 1 at the beginning of the simulation) and the mode fc* « \f\4>a is excited. Then at formation the wavelength of these 
perturbations is a fraction J?*/f?horizon = (ah* 1 )!! ~ 10~ 2 of the horizon size. Since inflation has ended, the comoving 
horizon ( aH )~ 1 is increasing, so aH is decreasing and the modes excited later will be an even smaller fraction of the horizon 


size. 




5 


Robertson-Walker (FRW) background are 


4> + 3 Htj) + A = 0 

(4) 

□X + 3 H~x + A x x 3 + 9 2 4> 2 X = 0 

(5) 


where H = a/a is the Hubble parameter, whose value is related to the total energy density p by the 
Friedmann equation 


H 2 


8itG 


( 6 ) 


In order to study the behavior of (j> and x that follows from the above, we will express x i n terms of modes 

Xk- 

x(t, i) = J d 3 k («(!)«-* + 4** W<***) ■ (7) 

The amplitude of 4> is still very large at the end of inflation, \ x x 2 g 2 (j) 2 , and Eq. (5) is approximately 
linear in We can then use the mode equation 

Xk + 3H\k + ^^2 + 9 2 Xk = 0 ( 8 ) 

to study the beginning of the reheating process. It will be convenient to introduce conformally rescaled 
fields (j> = a(j)/<j)Q, x = ax/<t> o- an< i time dp = dt/a and define a dimensionless time parameter and wave 
number 


t = VXcpor), k = k / (V\(j>o). 


(9) 


Following e.g. [11, 16], we study the field spectrum in terms of a comoving number density for the field x. 


n K 







( 10 ) 


and comoving energy density p K = ui K n K ,, where u K 


y^ 2 + ^efi = \J K 2 + q(t > 2 + 3(A X /A)x 2 . 


III. PREHEATING IN THIS MODEL 

We begin by briefly outlining some results from previous studies of preheating, beginning with the 
case A x = 0 (see [11] and references therein). We then use these to develop an approximate relation that 
quantifies the end of preheating and that will be useful in the gravitational wave calculation. After the end of 

7 Here we always use the Fourier Transform convention f(x) = (2n)~ 3 ^ 2 f d 3 k f(k) exp (ik ■ x). 
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inflation, 0 oscillates in its potential with period T ~ 7.416 (in terms of the dimensionless time parameter 
r) [11] while the modes Xn can be excited by the phenomenon of parametric resonance. This process is 
typically described in terms of the resonance parameter q = g 2 / A. In general, certain modes re will be 
excited as Xk oc exp(/i K r). The exponential growth factor ji r - will vary with re, giving rise to resonance 
“bands” characterized by some central re and width Are. We will consider the case of “broad resonance” 
where g > 1 (as compared with “narrow resonance” when q < 1). In this case the spectrum of resonantly 
excited modes takes the form of a broad peak whose location and width are approximately characterized by 


re*, Are ~ g 1 / 4 . 


( 11 ) 


For a particular value of q, the maximum growth exponent /i max = max{ p K } is [11] 

+ exp 


1 


/imax = - In ( \/l + exp -'KK 2 y/2jq 


-7rre 2 /v / 2g ) 




( 12 ) 


and the resonance is efficient when re 2 < \Jq/ (27r 2 ). Numerically we find that typical resonant momenta 
are re* ~ 1, so /x max ~ ( 3 / 27 t) exp(— iry/2/q) which is 0(1O _1 ) for the range of q we consider. Number 
density n x = f d 3 Kn XK increases in steps, twice per 0 oscillation - every time the inflaton passes through 
0 = 0 and x’s effective mass-squared m 2 = g 2 (jr goes to zero, a burst of \ particles are created. 

The exponential amplification of some \k derived from Eq. (8) is a solution for small % (approximately 
zero) and homogeneous 0, when the mode equation for is linear. As this process evolves, this will 
become a worse approximation and the problem will become fully nonlinear. Therefore, Eq. (8) is only 
useful for understanding the beginning of the reheating process, and in general it is the coupled equations 
of motion Eq. (4) and Eq. (5) that must be solved. 

These can be studied by lattice simulation, and we have used the C++ code LATTICEEASY [40] in or¬ 
der to simulate the evolution of these interacting scalar fields in an expanding universe. Fig. 1 shows results 
for q = 120. This is a useful example since [16] presents detailed results for preheating and gravitational 
wave production for q = 120 in the absence of a self-coupling. Fig. la shows the spatially-averaged energy 
density p x = x 2 + ^s(djX) 2 + \^\X A ) as a function of time. Fig. lc shows for A x = 0 the sum of the 
spatially-averaged energy densities p 0 = (^0 2 + ^(c0'0) 2 + qA0 4 ) and p x , as well as the energy density 
only in the interaction term, pi nt = (xfl 2 0 2 X 2 )- Fig- le shows the spectrum in x f° r the same choice of 
parameters. The spectrum is shown at several times, and the solid line corresponds to approximately the 
time when the exponential growth ends. 

We can understand how preheating progresses by observing that the transfer of energy between \ and 
0, and among different modes Xk and 0«-- occurs in the following distinct stages. First, oscillations of the 
homogeneous 0 excite modes of x centered around some re = re*, and the initially small inhomogeneities of 








L°H) / x djB L°<Px) / d.e 
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(a) (b) 



(c) 


(d) 



K K 


(e) 


(f) 


FIG. 1: Evolution of preheating for q = 120. (a) Energy density p x as a function of time for A x = 0. (b) 
Energy density p x as a function of time for A x = 10“ 8 . (c) Energy density of <i> and x, as well as energy 
density in the interaction term, for A x = 0. (d) Same as (c), but for A x = 10 -8 . The spatially averaged 
quantity q( jA x x 4 ) is also shown, (e) The spectrum in % at several times of interest, for A x = 0. The solid 
line corresponds to approximately the time when the exponential growth ends, (f) Same as (e), for 

A x = IQ" 8 . 
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X become large. There is some backreaction onto qb, whereby the /^-0 2 x” interaction term broadly excites 
modes <j> K up through ~ 2k,*, and inhomogeneities in <f> begin to grow. 

The second stage occurs once q 1//2 pi n t ~ p 0 + p x . This is a useful, approximate numerical result, that 
is essentially the same as Eq. 6 in [41]. Then x«* efficiently rescatters, i.e. interacts with other modes, 
and its exponential growth ends. The total energy in x continues to grow a bit until p x ~ p 0 . This is 
evident in Fig. la. Large field inhomogeneities break up and the spectrum broadens towards larger k. This 
broad spectrum where energy density becomes approximately evenly distributed among modes is evident in 
Fig. le. This figure indicates the spectrum at the time when the exponential growth ends with a solid curve. 
Spectra before this time are indicated by dashed curves, and spectra after this time are indicated by dotted 
curves. This stage is discussed and examples of field configurations are shown in [42]. Some work has also 
examined the final, so-called “turbulent thermalization” stage in detail [43, 44]. 

We now consider the case of nonzero A x . This has been studied to some extent in [34, 35, 45], and 
here we find results consistent with theirs. Fig. lb shows p x as a function of time. The resonance ends 
earlier in comparison with the A x = 0 situation of Fig. la. Fig. Id shows p, x + p x , p; nt and (jA x x 4 ) for 
A x = 10 -8 . Fig. If shows the spectrum in x at several times of interest, and the solid line again corresponds 
to approximately the time when the exponential growth ends. Here the end of this stage still corresponds to 
a large mixing between modes, but in this case it is the quartic self-interaction that is significant. 

In general, we find from numerical simulation that when A x becomes significantly larger than g 2 , the 
resonance terminates earlier than for the A x = 0 case, i.e. for any A x > A* ~ g 2 . In terms of energy 
transfer, when q 1 / 2 (|A x x 4 ) ~ Pint, the resonance ends. This is analogous to the condition we described 

for A x = 0, and will be useful. Depending on the size of A x , this may occur before or after the relation 

q 1 / 2 pint ~ Ptj, + p x becomes true. To summarize, the resonant stage of preheating ends by the following 
condition: 

{Pcfr + Px) ~ <? 1/2 Pint for A x < A*, (13) 

Pint ~ q 1/2 ^| A xX 4 ^> for A x > A*. (14) 

The powers of 1/2 are approximate - when comparing the size of the oscillating energy densities, as in 
Fig. lc and Fig. Id for example, there is some ambiguity in determining exactly what the value of the energy 
is when the resonance ends. We can estimate the value A* where the condition Eq. (14) becomes more 
important than Eq. (13) in terms of an energy argument. For small enough A x , we will have (jA x x 4 ) -C 
P<t> + P X . so fo° self-interaction will not play a role in ending the resonance. This will no longer be true once 



rsj 


p<t> + p x - 


(15) 
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This can be related to the value of x when the resonance ends by observing that, around this critical value 
A* where behavior transitions from Eq. (13) to Eq. (14), we will also have 

P4> + P X ~ 9 V2 ndXend^ (16) 

so that 

\\{xLd) ~ T^^endXend) ( 17 ) 

For ((0 2 X 2 >/(x 4 ))end ~ 0(1) this means that 

A* ~ g 2 . (18) 

This agrees with numerical results showing that the maximum energy density begins to decrease dramati¬ 
cally with increasing A x around this value. For example, q = 120 will give A* ~ 120 x 10 -13 ~ 10 _n . 
We check this by defining for each A x the quantity p™ ax as the time average over several oscillations once 
p x has stopped increasing with time. Fig. 2b shows that around A* ~ 10 ' 1 ’, p“ iax begins to decrease as 
A^ 1 . We now seek to quantify the effect that this has on gravitational wave production. 


IY. GRAVITATIONAL WAVE SPECTRUM 


The metric perturbation hij defined in Eq. (1) can be rescaled as h V] = ahij. Neglecting a term that 


goes as a"/a ~ ( aH ) 2 [16], the equation of motion is 


— V z hij = 167rGa 3 n,A 


(19) 


where G is Newton’s constant and 11 3 " 1 is the transverse traceless projection of the anisotropic stress: 


Ylij — a (Tij {p)Pij) ■ 


( 20 ) 


The second term in Eq. (20) will be neglected since g t j is the sum of a homogeneous, isotropic part whose 
transverse traceless projection is zero, and a perturbation that is higher order in G. The Fourier Transform 
of Eq. (19) is 

h”j(k) + k 2 hij(k ) = 16-7rGa 3 Uj^ (k) (21) 

We consider fI^ T to be a source acting continuously during the time interval r/o < >j < ?/j, solve 
Eq. (19) using Green’s functions, and use this solution to find the energy density of the tensor perturbation. 
As shown in [16], the result of this procedure is 


dpgw / \ Sk 


dink 


(V > Vf) = 


a 4 {v) 


( 22 ) 
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where Sk is defined by 


5fc = 4 ^ 3 [dn^2 


V 


rn f 


dr\ cos (kr/')a(r] , )T? T (rf, k ) 


’Vi 


+ 


rvf 


dr/ sm(kri , )a(r) , )T? T ( 7 /, k ) 


’m 


(23) 


where V is the volume of the box considered and J dfl is an integral over directions in k space. 8 Sk only 
depends on the dynamics occurring during gravitational wave generation, and the TT part of the energy- 
momentum tensor is defined in terms of projection operators by 



Pij(k) — dij 


h kt 


-Pij(k)Pi m {k) ) k) 


(24) 

(25) 


We obtain the spectrum of gravitational waves numerically using the LATTICEEASY code mentioned 
above, modified to in order to compute Eq. (22) as described above. We will give results in terms of ft gw = 
Pgw/ptotab at the “time of production” defined as approximately the time when energy in gravitational 
waves stops increasing noticeably. This is very well approximated by the value at the end of the simulation 
at r = 250, and denote with a subscript “p” quantities evaluated at this time. The relation between the 
results we give and their present values depends somewhat on the equation of state throughout reheating, 
but previous works have established that in A</> 4 preheating, the equation of state very rapidly becomes that 
of radiation, so that the energy density in gravitational waves will be [16] 

h 2 Q gw = (^) ^rad 

= (9.3 x 10" 6 ) (O gw ) p (26) 

where h 2 kl rM \ = 4.3 x 10" 5 , g*/go ~ 100. Similarly, frequencies today are related to comoving wave 
numbers at the time of preheating by 

f=( — ) 4 x 10 10 Hz ~ « x 10 7 Hz (27) 

\ap 1 PJ p 

where in the last step we have taken ( a 4 p) p ~ A0q (see e.g. Fig. la; we begin with p x « 0 and p^ « 
2 x jA(/>q and throughout the simulation the quantity « 4 (p d + p x ) ~ constant). 

Fig. 2a shows the spectrum obtained in the case q = 120, for A x = 0 and A x = 10 9 . The decrease 
in the energy produced in gravitational waves is evident from this, and Fig. 2b shows how this depends on 
A x , as a fraction of the peak energy density when A x = 0. The solid lines in Fig. 2c show lt gw for the cases 


Our physical results are independent of box size, as we use a numerical Fourier Transform that takes this into account. This is 
described in the LATTICEEASY documentation. 
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(C) 


(d) 


FIG. 2: Peak of gravitational wave energy density spectrum, defined in Eq. (22), as a fraction of total 
energy density at end of preheating stage, (a) Spectrum for q = 120 and two choices of self-coupling A x . 
(b) Amplitude of peak of GW spectrum, and final average value for p x after preheating ends, for q = 120 
and as a function of A x . These quantities are presented as fractions of their - value in the A x = 0 case. For 
comparison, dashed curves are also shown for the scaling behavior Eq. (28) and Eq. (29). (c) Amplitude of 
peak of GW spectrum, 0* w , as a function of A x for q = 12 and q = 1200, compared with Eq. (29). (d) 
Value of O gw as a function of the resonance parameter, q, for A x = 0 and A x = 10~ ', and the prediction 

Eq. (29) applied to the latter case. 


q = 12 and q = 1200. Evidently the effect of A x is to end the resonance early and suppress gravitational 
wave production. Once preheating ends, the additional contribution of inhomogeneities to the gravitational 
wave spectrum is negligible [16]. 

To estimate how this effect depends on the model parameters q and A x , we note that Q gw ~ (t£ t ) 2 ~ 
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(dix) 4 ■ The energy density Q gw is dominated by the most recently produced part of the spectrum before 
the resonance ends (this is particularly clear in Fig. 8 of [16]), so for the purposes of this estimate we will 
ask how the maximum amplitude of x depends on q and A x . We have seen that x grows until the condition 
Eq. (14), lqX(c/) 2 x 2 } ~ jA x (x 4 ), is satisfied. (Also, comparison of Fig. la with Fig. lb shows this since 
p x ~ (dix) 2 -) This suggests a parametric scaling 

Xend « Q X / X x = 9 2 / X x ( 28 ) 

Then the expectation that fi gw rsj (xLd) 2 becomes 

O gw oc ( 5 2 /A x ) 2 . (29) 

Our numerical results confirm this relation as shown in Fig. 2b and Fig. 2c. For A x < A*, the peak 
energy in gravitational waves decreases only very slightly with increasing A x , as the self-interaction term 
plays a small role in mixing modes and damping inhomogeneities. Once A x > A*, the energy density 
in gravitational waves scales in the manner given by Eq. (29). For A x ~ Iff 2 , we see that p x and Q gw 
no longer decrease significantly with increasing A x . This is simply because the unstable resonance never 
begins, and the quartic self-interaction can no longer dramatically decrease d gw by ending the resonance 
earlier. Fig. 2d shows how the value of O gw at the time of production depends on the resonance parameter, 
q, for both A x = 0 and A x = 10~'. In the latter case, we also show the prediction of the scaling relation 
Eq. (29). 


V. GENERALITY 

So far, we have examined results in the context of Ad 4 chaotic inflation, with the self-coupling A and the 
initial condition of the inflaton field identical to a previous work that thoroughly investigated the dynamics 
of gravitational wave production during preheating [16]. This allows the results of the previous sections to 
be directly compared with that work. However, observational data indicates that the Ad> 4 chaotic potential 
is not favored [37], so an important question is the generality of the results we have quoted above. In this 
section we will address this question in two ways, before pointing out interesting directions for future work. 
We will consider massive ( m 2 q i 2 ) inflation, another standard example in which preheating is studied, and 
we will also consider a range of initial conditions for cj) within both the Ad 4 and m 2 q i 2 cases. 

Specifically, this means that we will begin the numerical situation - corresponding to the end of infla¬ 
tion, with the inflaton’s energy about evenly split between kinetic and potential - with the inflaton field at 
various lower points on its potential than in the original case. Here, we are not primarily concerned with 
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FIG. 3: Peak of gravitational wave energy density spectrum, defined in Eq. (22), as a fraction of total 
energy density at end of preheating stage. Spectrum for q = 120 and two choices of self-coupling A x . The 
curves labeled 0(0) = 0o are identical to those shown in Fig. 2a, corresponding to the original initial 
condition for the inflaton field. The curves labeled 0(0) = 0o/lO correspond to starting the inflaton a 
factor of 10 lower on the potential, as described in the text. The magnitude of the gravitational wave 
spectrum is changed, but the effect of turning on A x is the same. 


representing a complete model of inflation, but rather are studying how preheating and gravitational wave 
production proceed within a potential that is quadratic or quartic about the minimum, without regard to 
the model’s behavior at higher (inflationary) field values. In this spirit, we also study the m 2 0 2 case with 
a few choices of m^. It is worth pointing out that not all inflationary models end with oscillations of the 
field responsible for inflation about its zero; see for example the Abelian Higgs and Higgs-dilaton models 
[20, 46]. 

For every situation we have tried, the same approximate scaling behavior of reduced gravitational 
wave production with increased self-interaction A x holds. In particular, we display some typical results 
in Fig. 3 and Fig. 4. For the case of A0 4 inflation, with q = 120 and 0(0) = 0o = 0.342 Mpi, we plot 
the gravitational wave spectrum in Fig. 3 and the scaling behavior with A x in Fig. 4a. These results were 
presented in Sec. IV, and they are provided again for direct comparison with alternative scenarios. We label 
this choice of parameters as 0 4 — I. We also show results for q = 120 and 0(0) = 0o/lO, referred to as 
0 4 — II, as well as q = 120 and 0(0) = 0o/lOO, referred to as 0 4 — III. 

Fig. 3 compares the gravitational wave spectra of the 0 4 — I and 0 4 — II parameter choices, for both 
A x = 0 and A x = 10 -9 . In both cases, evidently, there is a significant reduction in gravitational wave 
production that accompanies an increase in A x , despite the difference in overall amplitude of the spectrum. 
In Fig. 4a, we show how this reduction depends on A x for each of the parameter choices 0 4 — I, 0 4 — II, 
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(a) (b) 

FIG. 4: Peak of gravitational wave energy density spectrum, defined in Eq. (22), as a fraction of total 
energy density at end of preheating stage and normalized to the value when A x = 0. Here we show several 
typical examples where the initial condition and/or parameters of the model are varied, as described in the 
text. As in Sec. IV, there is some value A* above which the peak of the gravitational wave spectrum 
decreases as A~ 2 . Once A x is large enough, the preheating resonance never starts and there is no further 
suppression with increasing A x , an effect also seen in Sec. IV. (a) Varying initial conditions for Aci 4 
inflaton potential, (b) Varying initial conditions and mass parameter for m 2 (j) 2 inflaton potential. For ease 
of comparison, this result is given as a function of A x /A*, where A* = 10~ 6 , A 41 = 10~ 9 , A* 11 = 10 -7 . 

— III. We find the same scaling behavior as before: there is a A* above which gravitational wave 
production is suppressed by a factor of (g 2 /X x ) 2 . 

In the case of massive inflation, we replace Eq. (2) with the potential 

y = + \ x xX 4 + \g 2 4> 2 x 2 (30) 

i.e. the light field x has the same potential and interactions with the inflaton as it did previously, but the 
inflaton potential is quadratic rather than quartic. In this case we find that, as above, there is some A* such 
that for A x > A*, gravitational wave production tends to be suppressed by A x 2 . We again plot three typical 
examples. We refer to q = g 2 (f)(0) 2 /4m^ = 60, (f)(0) = 0.1 Mp\, = 10 -6 Afpi as 4> 2 — I. We refer to 
q = 15, 0(0) = 0.01, rriff, = 10 -9 as cf) 2 — II. We refer to q = 15, (f)(0) = 0.001, = 10 -9 as cf> 2 — III. 

Fig. 4b shows how the gravitational wave spectra in these cases scale with A x . For ease of comparison 
with the scaling relation A x 2 we plot the results as a function of A x /A*, where A* = 10~ 6 , 10~ 9 , 10 - ' 
for cf) 2 — I, (f> 2 — II, <f) 2 — III respectively. As before, A^ 2 fits well (until A x becomes large enough that 
preheating no longer starts, so that increasing A x won’t further decrease the gravitational wave production). 











15 


The numerical computations involved make it impractical to check here every imaginable situation of 
interest to verify this relation. We have shown that gravitational wave production from preheating in poten¬ 
tials with minimum at zero can be extremely sensitive to the value of the light field’s self-coupling term, 
and that result is not exclusive to one particular model or choice of parameters. Therefore, an important 
goal of future work will be to fully characterize this effect in other realistic models, and better understand 
the implications for observability. 


VI. DISCUSSION AND CONCLUSIONS 

In this work we have studied the effect of a nonzero self-interaction on gravitational wave production 
during preheating of a scalar field. Previous work has considered the dynamics of preheating for a light, self¬ 
interacting scalar, as well as gravitational wave production by preheating of a non-self-interacting scalar. 
This work is an extension of these results, and in particular shows that the spectrum of gravitational waves 
that survive until today is very sensitive to the light scalar’s self-interaction. Our main result within the A</> 4 
model is that for self-coupling A x > g 2 , the preheating resonance is terminated early, and the gravitational 
wave spectrum is significantly reduced: 

/ 2 \ 2 

ttgw ~ ^w x=0) for A x > g 2 . (31) 

We have also begun to address the question of generality of this result, as discussed in Sec. V. For 
various choices of the inflaton’s initial condition in the A<5 4 model, we have seen that Eq. (31) holds. 
Additionally, for an m 2 <j) 2 inflationary potential, the result that the gravitational waves are suppressed as A” 2 
is shown, for several parameter choices. While this suggests generality to inflation models with potentials 
quadratic or quartic about a minimum at zero, an important question for future work is to study the effect of 
the light field’s interactions in other preheating models that have been shown to predict gravitational waves. 
As our work shows, predictions that neglect such interactions - even if they are extremely small - may not 
necessarily be accurate. 

It is easy to imagine that in a realistic preheating scenario, decay products will have their own self¬ 
interactions or further interactions with other fields, that will end the resonance early. Recently, another 
paper studied the effect of interactions of x with further light degrees of freedom, as well as self-interactions 
in the context of a curvaton decaying to Higgs [32], Although the model is not identical to ours, it also found 
that self-interactions can be important in terminating the resonance early. Furthermore, they found that 
interaction with the additional light scalars, as characterized by the contribution to a thermal term, has the 
ability to significantly affect the resonance and either end it early or prevent it from occurring at all. They did 
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not consider gravitational wave production, but following the argument given here in Sec. IV it is reasonable 
to expect that this early termination of the resonance can further reduce any production of gravitational 
waves. Analyses of other scenarios have shown that preheating can be sensitive to nonlinear interaction 
terms of decay products [23], or other nonperturbative effects motivated by new physics above the TeV scale 
[47-51]. Another interesting goal for future work would be to incorporate the effects of interactions such 
as those studied in this paper into a more general framework for obtaining order-of-magnitude estimates 
of gravitational wave production, as in [33]. Although current constraints on MHz gravitational wave 
backgrounds are not sensitive to these processes [29], this could be very useful in evaluating the potential 
for observability in future experiments. 

One interesting possibility is that reheating occurred through an inflaton-to-Higgs coupling, since the 
Higgs is a natural candidate to couple to beyond-Standard Model fields [52-55]. The running of the Higgs 
self-coupling is sensitive to any new physics that comes in at high energies, but it has not been directly 
measured and will be difficult to measure at the LHC. One might hope that since Ah runs from 0.13 at 
the weak scale to zero around 10 9 10 * GeV in the Standard Model [56], the condition Ajj -C 1 could be 
satisfied. This would avoid enormous damping of the preheating resonance, and thereby provide a possible 
cosmological probe of Ah- The self-coupling remains O (l0 _1 ) up to ~ 10 8 GeV, though, which suggests 
that there will not be significant (or any) preheating resonance. However, above this scale the self-coupling 
decreases and the effective potential reaches a maximum (in the Standard Model - small changes in input 
parameters or new physics beyond the Standard Model can significantly affect this; see e.g. [57-60]). 

The condition Eq. (14) suggests that a more relevant condition than the self-coupling may be the mag¬ 
nitude of the Higgs potential. The configuration of x at the end of inflation (initial configuration for this 
problem) is certainly sensitive to the potential at large field values, as it corresponds to approximately 
Arms ~ 10 12 GeV Tfinf. 9 If one takes Eq. (14) to apply as the condition for whether parametric reso¬ 
nance does or does not occur, then the result could be a resonance pushing Higgs oscillations toward the 
vacuum instability region. 10 New physics that prevents Ah from becoming negative would likely be more 
than sufficient to prevent a resonance from occuring. These rough estimates also ignore the possibilities of 
a different running of Ah from the new inflaton coupling, as well as thermal effects. We leave the resolution 
of these questions to future work. 


9 The behavior of the Higgs after inflation, when there is no coupling to the inflaton, is discussed in [61]. 

10 There has been much work on the Higgs and vacuum stability, including discussion of the reheat temperature; see e.g. [62-64] 

and references therein. 
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